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ABSTRACT 

Ionizing sources embedded in the neutral intergalactic medium (IGM) before cosmological reionization 
generate discrete HII regions. We show that a sufficiently bright quasar (for example, one tenth as 
luminous as that recently discovered by Fan et al.) can ionize a large volume, allowing the transmission 
of a substantial fraction of the flux of its intrinsic Lya emission line on both sides of the Lya wavelength. 
The observed line profile is richly informative. We show that a sufficiently accurate, high spectral 
resolution (R « I0 4 ) measurement of the line profile of a bright quasar is feasible using the Next 
Generation Space Telescope (NGST) as well as large ground-based telescopes. Such a measurement has 
two potentially important applications. First, the red side of the Lya emission line provides a way to 
measure the quasar lifetime. Second, the blue side provides a direct measure of the density fluctuations 
of the intergalactic medium at the quasar redshift. The estimate of the absorption of the red side is 
limited by the accuracy to which the intrinsic profile of the Lya emission line is known. The blue side, 
however, does not sensitively depend on the intrinsic profile, because the former is much narrower than 
the latter. 

Subject headings: Line: profiles - quasars: absorption lines - quasars: emission lines - techniques: 
spectroscopic - ultraviolet: general 



1. INTRODUCTION 

The epoch of cosmological reionization, and the nature 
of the ionizing sources bear fundamental cosmological im- 
portance, but are yet unknown. In contrast, we have a 
better understanding of how the reionization might have 
happened: individual HII regions, whether around galax- 
ies or quasars (or sources of other nature), start out as 
isolated islands, eventually overlap to fill the universe and 
complete the reionization process (Aaron & Wingert 1972; 
Shapiro & Giroux 1987; Gnedin & Ostriker 1997; Haiman 
& Loeb 1998; Miralda-Escude et al. 1999). In this Letter 
we suggest that it may be possible to directly probe the 
structure of an HII region around a bright quasar before 
the reionization is complete. Although we do not necessar- 
ily have to identify the primary reionization sources with 
quasars, the existence of bright quasars at high redshift is 
supported by the recent discovery of the very bright quasar 
at z — 5.8 in the Sloan Digital Sky Survey (SDSS; Fan et 
al. 2000; F00 hereafter). 

To study the size and structure of an individual 
Stromgren sphere, we propose to utilize the observed pro- 
file of the quasar's Lya emission line. When a source is 
embedded in a neutral IGM, its emission line suffers from 
strong absorption by the damping wing of the resonance 
Lya absorption of the intervening neutral IGM (Miralda- 
Escude 1998). As a result, the Lya photons are absorbed 
and re-emitted multiple times in the IGM, resulting in an 
extended, low-surface brightness line, which is difficult to 
observe (Loeb & Rybicki 1999). In contrast, if the source 
is surrounded by a large HII region, the Lya photons can 
escape absorption, provided they redshift out of resonance 
before they reach the boundary of the HII region. In this 
Letter, we show that a sufficiently bright quasar placed at 
z > 6 can produce a sufficiently large Stromgren sphere 



to render its Lya emission line detectable, using the Next 
Generation Space Telescope (NGST), as well as high spec- 
tral resolution ground based telescopes such as Keck and 
the Hobby-Eberly (HET) telescopes. Throughout this Let- 
ter, we adopt a flat ACDM background cosmology dom- 
inated by a cosmological constant and cold dark matter, 
with density parameters (f2A,^o,^t>) = (0.7,0.26,0.04), 
and Hubble constant Ho = 70 km s . 

2. STROMGREN SPHERES AROUND QUASARS 

According to the classical definition in stellar astronomy, 
a Stromgren sphere is an HII region around an individual, 
ionizing-photon-cmitting O/B star, within which the rate 
of recombinations exactly balances the emission rate N p h 
of ionizing photons from the star. The radius R s of the 
Stromgren sphere is given by 
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where cub is the hydrogen recombination coefficient at 
T = 10 4 K, and (n 2 H ) is the mean squared hydrogen density 
within R s . 

An analogous situation arises on much larger scales, for 
cosmological ionizing sources such as galaxies and quasars 
(Shapiro & Giroux 1987). There are, however, important 
differences between the stellar and cosmological case. The 
most important one for our purposes is that the hydrogen 
recombination time in the IGM is longer than the Hubble 
time at the relevant epochs (z <, 10) and the quasar life- 
time. As a result, the size R s of the equilibrium Stromgren 
sphere is large, but is never reached in practice before the 
ionizing source turns off. The radius Rt Q of the ionized 



bubble Fellow 



2 



region around a steady source of age tQ is then simply 
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where (tih) is the mean hydrogen density within Rt Q - A 
second difference is that in the stellar case, the Stromgren 
sphere is static if the luminosity of the source is constant, 
whereas in the quasar case the sphere around a long-lived 
quasar can evolve due to the evolution of density around 
a quasar (both mean density and clumping factor) as well 
as the accumulation of ionizing photons (equation 2 and 
equation 4 below), even if the source luminosity remains 
constant. To fully take these effects into account, one 
needs to solve the equation of motion for the ionization 
front, Ri, with the following equation 



dt w ' 47r(n ff ) 
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where H(z) is the Hubble constant at z, and Chii = 
(n 2 H ) I '(nn) 2 is the mean clumping factor of ionized gas 
within Ri. The three terms on the right side in equa- 
tion [| account for the Hubble expansion, the ionizations 
by newly produced photons, and recombinations (Shapiro 
& Giroux 1987; Haiman & Loeb 1997), respectively. Note 
that equation |3| predicts an initial phase of superluminal 
expansion for the I-front, i.e. when Ri(t) > ct. Ignoring 
any clumping or recombinations, this phase would last for 
t » 10 7 (iV p/t /10 57 s- 1 ) 1 / 2 [(l + z)/8}- 3 / 2 years (cf. eq. §); 
a potentially significant effect, since the quasar found by 
F00 in the SDSS survey implies an intrinsic ionizing emis- 



sivity of N ph « 2 x 10 



57 



Nevertheless, the Lya pho- 



tons emitted by the quasar lag behind the I-front during 
the phase of relativistic expansion, and will only see the 
opacity once the front slows down to be subluminal. As a 



result, for our purposes, the a 
is always given by equation B 
In practice, we find that t 



^parent size of the HII region 

le apparent size of a quasar 
Stromgren sphere is primarily determined by the total 
number of ionizing photons emitted. Assuming that the 
quasar luminosity N p h is constant, the solution of equa- 
tion H is accurately described by R tQ given in eq. |[ Gas 
clumping is unlikely to significantly change this answer. 
As a specific example, with the ionizing photon emissiv- 
ity of the SDSS quasar, and for Chii = 100, we find that 
Ri would be reduced by 10% at « 10 7 yr, and by 50% at 
« 10 s yr. Cosmological simulations indicate that the value 
of Chii is likely to be significantly below 100 (Gnedin & 
Ostriker 1997). 

3. LYa TRANSFER ACROSS A STROMGREN SPHERE 

The main purpose of this Letter is to assess the effect of 
the quasar's HII region on the observed profile of its Lya 
line. Consider a source located at z s <; 6 before reioniza- 
tion is complete, whose age is tQ. The source is assumed 
to be surrounded by a spherical HII region of size Ri(tQ), 
but is otherwise embedded in a neutral IGM. The source 
is also assumed to emit a Lya emission line, which is sub- 
sequently reprocessed by the opacity of the intervening 

2 We assume that the IGM temperature is coupled to the cosmic 



neutral IGM, as well as of the residual neutral hydrogen 
within the HII region itself. 

The optical depth between the source and the observer 
at z = 0, at the observed wavelength A Q b s = A s (l + z s ), is 
given by 
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where cdt/dz is the line element in the assumed ACDM 
cosmology, nn is the neutral hydrogen density, and a a is 
the Lya absorption cross-section. The reionization red- 
shift is assumed conservatively to be z r = 6, but the pre- 
cise value (as well as the details of the reionization process) 
is irrelevant here, since we are only sensitive to the opti- 
cal depth near the source redshift z s , where the IGM is 
assumed to be still fully neutral. It is useful to divide 
the range of integration in equation |] into two parts, to 
reflect the contributions to the optical depth from within 
(zi < z < z s ) and outside (z r < z < Zi) the HII region, 
where Zi is the redshift somewhat below z s , corresponding 
to the boundary of the HII region (zi w z s — R4 [t] /Rh [z s ] , 
where Rh[z s ] is the size of the cosmological horizon at z s ). 

The Lya absorption cross-section is assumed to be de- 
scribed by a numerically computed Voigt profile (see eq. 6 
of Press & Rybicki 1993). The gas temperature for this 
profile is taken to be T = 10 4 K within the HII region, and 
T = 1.2[(l + z)/8]~ 2 K outside it. 2 Outside the HII region, 
the neutral hydrogen density is assumed to follow the IGM 
density, 8.5 x 10" 5 [(1 + z)/8} 3 cm -3 . Inside the HII re- 
gion, the neutral hydrogen density depends on the ionized 
fraction. Assuming photoionization equilibrium (justified 
for the source lifetimes of > 10 6 yr considered here), and 
that the gas is highly ionized, the neutral hydrogen frac- 
tion x ~ riHi/nn C 1 at the quasar-centric radius r is 
given by 
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Here L v is the intrinsic luminosity of the source in 
erg s^ 1 Hz -1 , vjj and a v are the ionization threshold and 
cross section of hydrogen, with the luminosity weighted 
cross section being 
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and equal to 2.5 x 10~ 18 cm 2 for the Elvis et al. (1994) 
template spectrum (the same value would be obtained for 
a power law spectrum with a slope of v" 1 - 8 ). The resulting 
neutral hydrogen fraction is 
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This equation applies within the optically thin central re- 
gion of the Stromgren sphere where a; C 1, and shows that 
for bright sources of interest here, this region is highly 
ionized. Near the boundary of the HII region, the neu- 
tral fraction rapidly increases to unity, where the x <C 1 

microwave background until z = 150, and cools adiabatically thereafter. 



condition breaks down and equation ^ becomes invalid. 
However, this does not affect the results shown below for 
the following reasons. First, the flux transmission short- 
ward of the Lya wavelength that we are concerned with 
occurs in regions with optical depth r<7ora;<5x 1CP 5 . 
Second, the optical depth longward of the Lya wavelength 
is roughly 0.5a;(AA/100A ) _1 for the cosmological param- 
eters adopted here (Miralda-Escude 1998), and is dom- 
inated by neutral (x = 1) IGM outside the Stromgren 
sphere. Therefore, the boundary region of the Stromgren 
sphere does not affect the results in a sensitive way. We 
assume that the jump from x(Ri) to unity takes place as 
a step function at the radius Ri(t,Q). 



it is quite insensitive to £q (equation [?]) . 
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Fig. 1. — The total line optical depth (solid curves) towards 
a quasar at z a = 7, as a function of observed wavelength. Also 
shown are the separate contributions from the neutral IGM (dot- 
ted curves) and from the residual neutral hydrogen within the HII 
region (dashed curves). In the left panel, we assume an ionizing 
photon emissivity of 2 X 10 57 s — 1 , a fixed quasar lifetime of 10 7 yr, 
and demonstrate the effect of the clumping factor (Chii = 3, 10, 30). 
Clumping effects primarily the blue side of the observed Lya line. 
In the right panel, we assume an ionizing photon emissivity of 
2 X 10 5e s _1 , a fixed clumping factor of Chii = 10, and demonstrate 
the effect of the quasar lifetime (10 6 , 10 7 , 10 s yr). The lifetime ef- 
fects primarily the red side of the observed Lya line. 



4. RESULTS AND DISCUSSION 

Figure [l] shows the optical depth around the Lya line 
for a quasar at z s = 7. The left panel shows that, for 
a bright quasar like that of F00, the optical depth long- 
ward of the Lya wavelength is nearly independent of the 
clumping factor, simply because the optical depth in that 
region is dominated by the damping wing of the neu- 
tral IGM (dotted curves in the left panel) outside the 
Stromgren sphere, whose radius is insensitive to Chii as 
long as Chii < 100 (see §2). The situation is dramatically 
different for a fainter quasar as shown in the right panel 
of Figure [IJ because the HII region is much smaller, the 
red damping wing produced by the neutral IGM outside 
the Stromgren sphere start to reflect the size of the sphere 
hence the age of the quasar. The dependence of the red 
damping wing on the clumping factor is very weak, be- 
cause the contribution from the residual neutral hydrogen 
within the Stromgren sphere is small. 

On the other hand, the optical depth shortward of the 
Lya wavelength within the stretch of the Stromgren sphere 
is always dominated by the residual neutral hydrogen in- 
side the Stromgren sphere (dashed curves), and therefore 
is proportional to Chii (equation fj]), in sharp contrast with 
the situation at the longward of the Lya wavelength. How- 
ever, the optical depth shortward of the Lya wavelength 
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Fig. 2. — The transmitted line profile for the Lya line of a quasar 
with a central black hole mass of either 2 X 10 9 Mq (similar to what 
is inferred for the z = 5.8 SDSS quasar found by F00, upper panels), 
or 2 X 10 s Mq (lower panels) . A source redshift of either z s = 10 (left 
panels) or z s = 7 (right panels) is assumed. The upper dashed curve 
shows the line profile without any absorption. The set of three solid 
curves shows the transmission assuming a quasar lifetime of 10 7 yr, 
and using increasing clumping factors, corresponding (top to bot- 
tom) to Chii = 1, 3, 10 (z = 10) and Chii = 3, 10, 30 (z = 7). The 
dotted lines show the same results, except for a quasar lifetime of 
10 8 yr. The lower dashed curve assumes a quasar lifetime of 10 6 yr, 
and Chii = 1 at z = 10 and Chii = 1 at z = 7. Also shown is the 
level of noise expected in a 10 5 s integration with an 8m telescope, 
such as NGST. 

Let us now translate the theoretical illustration in Fig- 
ure [l] to observable spectra of quasars around the Lya 
emission line shown in Figure p[ A logarithmic ordinate is 
used to better show the flux over many orders of magni- 
tude down to the expected noise level (although the spec- 
tral line may appear misleadingly flat). We adopt the 
black hole (BH) mass and luminosity inferred from the 
quasar at z = 5.8 observed by F00 and place it either at 
z s = 10 or z s = 7 in the top two panels; another quasar 
with one tenth of this luminosity is used for the two bot- 
tom panels. In all cases, the profile of the Lya emission 
line is assumed to be Gaussian with a half- width of 1 , 500 
km/s, and with a central flux that is twice the continuum. 
The random noise shown at the bottom of each panel in 
Figure || is estimated assuming a 10 5 second integration 
with the NGST (see bttp: / / augusta.stsci.edu for details 
on the expected backgrounds) with a spectral resolution 
of R=10,000. Obviously, it is advantageous to use high 
resolution spectroscopy to maximize the number of pixels 
that probe the structure of the HII region on the blue side 
of the emission line (especially for faint and/or very young 
quasars). At spectral resolutions this high, the dominant 
absolute (not relative) noise is either the Poisson photon 
noise (for bright sources with N p h 10 56 s _1 ) or the de- 
tector noise (for fainter sources), rather than the zodiacal 
light or sky background (see expressions for noise in Gillctt 
& Mountain 1998). A ground based telescope such Keck 
or HET with sufficient spectral resolution should therefore 
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work just as well as NGST. 

A rich amount of information is contained in the profile 
of the line. We see that the red side of the Lya emis- 
sion line is largely transmitted even for the fainter quasar 
adopted at the bottom panels, consistent with Figure 1 
showing an optical depth of ~ 0.1 there (except near the 
Lya wavelength for the case with tQ = 10 6 yr, where the 
radius of the Stromgren sphere is sufficiently small that the 
damping wing of the neutral IGM contributes an optical 
depth significantly above unity.) For the brighter quasar 
shown in the top panels the transmission on the red side 
is insensitive to either tQ or Chii, whereas for the fainter 
quasar shown in the bottom panels the transmission on 
the red side is sensitive to <q but not to Chii- It is there- 
fore possible to measure the lifetime of a relatively faint 
quasar (perhaps the majority of the quasars belong to this 
category) before reionization is complete. The primary 
uncertain factor is the intrinsic Lya emission profile. 

The blue side of the Lya emission line is not completely 
absorbed. On the contrary, we see in the top panels of 
Figure || that for a quasar observed by F00 with tQ > 10 7 
yr there is a large range A A w 50 — 150 A where the flux 
is at least ten times above the noise (S/N=10). Even if 
the quasar is ten times fainter, one should still be able to 
detect of order several tens of spectral pixels if tQ > 10 6 
yr (assuming R = 10 4 ). 

More interestingly, the extent of the flux transmission 
T v on the blue side is a strong function of Chii but a 
very weak function of tQ . The idealized case in Figure |^ 
shows the mean flux processed through a relatively gen- 
tly fluctuating medium, with no saturated lines [so that 
T v oc exp(— /3i/Chii), where (i\v is a known constant of 
order 0.1 — 1.0]. An actual observed flux distribution on 
the blue side may contain numerous absorption features 
due to density fluctuations on small scales. In the limit 
when clumping is caused by highly overdense regions with 
a small volume filling factor, the absorption lines saturate, 
and the mean transmission could increase, rather than de- 
crease, for increasing values of Chii (since a larger fraction 
of the line of sight probes underdense voids). In order to 
adequately model the absorption profile and to make use 
of the additional available information, numerical simula- 
tions are needed, which we will explore in a subsequent 
paper. In either case, the blue side of the Lya line will 
provide a potentially powerful tool to directly measure the 
small scale power at the quasar redshift. 

Since the intrinsic Lya emission of a quasar is expected 
to be significantly broader than the width of the blue side 
of the absorbed spectrum, the latter should not sensitively 
depend on the exact intrinsic width of the unabsorbed Lya 
emission line. Since the red side of the profile traces the 
intrinsic profile well for a bright quasar, it is possible that 
fitting the entire observed profile significantly removes the 
uncertainty of the unknown intrinsic profile, as long as it 
is reasonably symmetric and sufficiently broad. 

Hogan, Anderson, & Rugers (1997) and Anderson et al. 
(1999) have made the point that the quasar lifetime and 
ionizing background from the proximity effect for a HcIII 
region around the quasar (at lower redshift) is related and 
showed that a large Helll region can be created by a bright 
quasar with a reasonable lifetime. They indicated that for 
a bright quasar of age tQ ~ 10 7-8 yr a large Helll region 



of size ~ 20A can be generated. 

Finally, we point out that the existence of a useful HII 
region around a quasar requires only that the reionization 
is not complete. When individual HII regions overlap, the 
universe is not yet completely reionized, due to the travel 
delay of ionizing photons from other sources. In other 
words, the emission profile that we show above should 
exist even at times significantly later than the overlap- 
ping epoch, but before the average Gunn-Peterson optical 
depth of the IGM is reduced to below unity (note that the 
latter can take a significant time, see, e.g. Gnedin 2000). 
This could potentially significantly reduce the redshift of 
quasars that should show useful transmission profiles. 

5. CONCLUSIONS 

A quasar one tenth as luminous as that recently re- 
ported by F00, if placed in the neutral IGM at z s > 7 
before cosmological reionization, generates a large HII is- 
land (a Stromgren sphere). We show that for such a quasar 
a substantial fraction of the flux on both the blue and 
red sides of the Lya emission wavelength can be transmit- 
ted. For the range of wavelengths around the Lya wave- 
length where significant transmission of flux occurs, the 
red damping wing of the intervening neutral IGM outside 
the Stromgren sphere dominates the optical depth long- 
ward of the Lya wavelength, whereas the residual neu- 
tral hydrogen inside the (highly ionized) Stromgren sphere 
dominates the optical depth shortward of the Lya wave- 
length. While the flux longward of the Lya wavelength is 
not substantially suppressed to a large offset in wavelength 
(SX > 100A), the blue side of the Lya emission line is more 
absorbed than the red side and displays a distinct profile 
that drops sharply with decreasing wavelength. Neverthe- 
less, the flux on the blue side is clearly detectable even for a 
quasar 10 times fainter than that found by F00. We show 
that NGST as well as large ground-based, high spectral 
resolution telescopes with a spectral resolution of R = 10 4 
should be able to provide such a detection. 

Two applications of such a partially absorbed Lya emis- 
sion line are suggested. First, for a relatively faint quasar 
(for example, one tenth as luminous as that found by F00), 
the red side of the Lya line profile depends mostly on the 
age of the quasar and therefore can be used to infer it. The 
main uncertainty is the unknown intrinsic profile. Second, 
for bright quasars, the blue side of the Lya line profile 
is also detectable, and depends only very weakly on the 
lifetime of a quasar but very strongly on the density fluc- 
tuations the gas around the quasar. Therefore, one can 
use it to directly measure the fluctuation of gas on small 
scales at high redshift. This may provide a unique oppor- 
tunity to probe the small-scale power at very high redshift, 
where competing cosmological models differ most. 
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